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Abstract. - Magnetic nanowires are a good platform to study fundamental processes in 
Magnetism and have many attractive applications in recording such as perpendicular storage 
and in spintronics such as non-volatile magnetic memory devices (MRAM) and magnetic logic 
devices. In this work, nanowires are used to study magnetization reversal processes through a 
novel geometric approach. Reversal modes imprint a definite signature on a parametric curve 
representing the locus of the critical switching field. We show how the different modes affect the 
geometry of this curve depending on the nature of the anisotropy (uniaxial or cubic anisotropy) , 
demagnetization and exchange effects. The samples we use are electrochemically grown Nickel 
and Cobalt nanowires. 



Introduction. - A surge of renewed interest in ferromagnetic nanowires has occurred very 
recently triggered by their interesting properties with respect to spintronic devices and more 
specifically to non- volatile memory (MRAM) and logic devices [1-3]. 

From a fundamental point of view, they represent a quasi-one dimensional system dis- 
playing magnetic properties in sharp contrast to the Mermin- Wagner [4] theorem forbidding 
(Hcisenberg-type) Magnetism in systems of dimension (D < 2) with short-range interactions. 

They are simpler than nanotubes since their physical properties do not depend on chirality 
and they can be grown in a variety of methods: Molecular Beam Epitaxy, Electrochemical 
methods (Template synthesis, Anodic Alumina filters), Chemical solution techniques (Self- 
assembly, Sol- Gel, emulsions...) etc... 

Template synthesis is an electrochemical method used in the manufacturing of nanostruc- 
tured materials, in particular, nanowires can be obtained by filling a porous polycarbonate 
membrane, wich contains a large number of cylindical holes that are track-etched and pos- 
sessing a narrow size distribution. 

Individual nanowires as well as arrays of nanowires are important. Characterization and 
understanding of the magnetic properties of nanowire arrays are challenging since we have 
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dipolar interactions between wires in addition to anisotropy, magnetostatic and exchange 
effects within individual wires. 

Many questions remain open regarding the detailed mechanisms responsible for the mag- 
netization reversal. The intrinsic properties of nanowire arrays are directly related to the 
properties of the nanoporous template such as the relative pore orientations in the assembly, 
pore size and its distribution, as well as interpore distance. 

Two reversal modes have been suggested as being important: curling and the coherent 
rotation depending on the value of the radius a with respect a cutoff value also called " coher- 
ence" radius R co h- 

For wire radius a > R co h, the reversal occurs by curling. For a < R co h, the coherent 
rotation is expected as predicted by the Stoner-Wohlfarth model [5] . 

The remanent magnetization distribution within a nanowire depends on its radius as well 
as on the balance between different contributions (anisotropy, magnetostatic/shape and ex- 
change) to the magnetic energy. Exchange energy dominates at small wire radius and favors 
a uniform magnetization distribution or non-uniform flower, S, Landau, leaf... states. On 
the other hand, this almost collinear spin alignment leads to large demagnetizing fields due 
to magnetic surface charges and correspondingly increase of the magnetostatic energy which 
increases gradually with radius. When the radius increases beyond the single domain limit 
(a > R s d), the competition between these two energies and anisotropy energy leads to non- 
uniform magnetic states such as multi-domain states (if the anisotropy constant is greater 
than the magnetostatic energy). It might lead to magnetic vortices, if the anisotropy energy 
is small in comparison with magnetostatic energy. 

In this work, we set out to study the geometric signature of reversal modes in nanowires 
depending on the different magnetic properties inherent to the wire belonging to any of the 
anisotropy, magnetostatic and exchange classes. This paper is organised as follows: In section 
II, the geometric signature is discussed in the case of an infinite cylinder with a uniaxial 
anisotropy. In section III, the resulting parametric curve is considered for a general ellipsoid 
with arbitrary uniaxial anisotropy and finally section IV bears our conclusions. 

Geometric signature of an infinite cylinder. The geometric signature of a reversal 
mode is a parametric curve derived from the nanowire magnetization configuration that can 
be determined from the Brown equation by minimizing the total free energy [6] . 

In bulk ferromagnetic materials, the energy of the system can be minimized by forming 
multiple magnetic domains within which the magnetic moments are aligned. However, there 
is a critical size R s d below which a particle remains in a single-domain state during switching; 
it is given by the implicit form: 

derived by Frei et al. [7] after approximating the nanowire by an ellipsoid of revolution 
with major axis c and minor axis, a = b / c and comparing the exchange energy averaged 
over the ellipsoid volume to the average magnetostatic energy (see Fig. 1). 

For Nickel, the following parameters might be used: M s = 485 emu/cm 3 , ao = 0.249 
nm [8]. 

The curling reversal mode is the dominant magnetization reversal process in magnetic 
nanowires. The magnetization curling mode was defined by Frei et al. [7] and after it has 
used for different structure to investigate the magnetic switching of films, spherical particles, 
prolate ellipsoid and cylinders [9] . 
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Fig. 1 - Cylindrical nanowire approximated by an ellipsoid of revolution of radius a = b and long axis 
c. The nanowire is considered as a single domain with uniaxial anisotropy and in presence of applied 
field H. The long axis of the ellipsoid taken along the z direction is also the anisotropy axis. The 
(u, w, y) frame is such that the equilibrium magnetization Mo direction defines the u axis. Mo lies 
in the (y, z) plane as well as the (u, w) basis obtained from it through rotation by the equilibrium 
angle 8 (see Appendix). 

The coherence radius R co h separating uniform rotation and curling is given by: 



A is the exchange stiffness (^4 = 1.5 x 10~ 6 erg/cm for bulk Nickel and Cobalt), M s is the 
saturation magnetization and N a is the demagnetization coefficient along the nanowire radius 
(see Fig. 1). The parameter q, is the smallest positive zero of the first kind Bessel function 
J\ (x) derivative dJ ^f* [10] (i.e. smallest positive maximum) and the smallest of the derivative 
of all ellipsoidal harmonics. It has been evaluated by Aharoni [11] to a 10~ 7 accuracy. He 
also found an accurate interpolation formula for any aspect ratio m — c/a: 



with coefficients : b = 1.84120, bi = 0.48694, b 2 = -0.11381, b 3 = -0.50149, 6 4 = 
0.54072, 6 5 = -0.17200, that are accurate to 10~ 6 . 

For particle size larger than R co h yet smaller than R sc i, magnetization reversal proceeds 
through curling. In this mode, magnetization switching is an abrupt process, and the switching 
field is very close to the nucleation field for all angles. 

In the case of an Infinite cylinder with a uniaxial anisotropy along the axis K, two distinct 
calculations were made independently by Chung and Muller [12] and Ishii [13] 20 years apart. 

The curling equations (see Appendix) to be solved are: 
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It is remarkable that the constant k 
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that appears in a = is a rational 

approximation to the first zero of the Bessel function derivative dJ ^ x ^ ■ This originates from 
the way Ishii [13] used the Bessel function in the average energy calculation (see Appendix 
and next section). The parameter S is the reduced radius defined as S = a/£ ex , with £ ex the 
exchange length [14] considered as an intrinsic length scale of the nanowire (e.g. £ ex —20.6 
nm for Nickel). £ = in coherent rotation (Stoner-Wohlfarth [5] limit), 5 = 1 in buckling 
whereas S > 1 in curling with the corresponding parametric curves displayed in Fig. 2. 
The switching field components h x = H 2 ^^ , h z = H 2 c °m^ are then obtained as: 



h x = sin 9 



1 -a + 
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The above equations do not produce the coherent rotation limit as S 
one writes: 



(5) 



in which case, 



^ = -(^ + l)cos 3 * (6) 

There is a discrepancy in the results obtained by Ishii [13] and Chung and Muller [12] 
whose a and anisotropy A factors are off by a factor of ^. 

Geometric signature of a finite ellipsoid with arbitrary uniaxial anisotropy. Let us 
consider the general case of an ellipsoid with fourth-order uniaxial anisotropy characterized 
by two constants Ai and Following Aharoni [9], Wegrowe et al. [17] performed a detailed 
study of Template synthesised Nickel nanowires and concluded that the anisotropy of their 
nanowires is positive and large (recall that for bulk Nickel, the lowest order cubic anisotropy 
constants [15] are Ai = —4.5 x 10 4 erg/cm 3 and K 2 = 2.3 x 10 4 erg/cm 3 at room temperature) 
arguing that a large strain acting on the nanowires may induce such a large anisotropy change. 

The curling equations to be solved in presence of fourth-order uniaxial anisotropy are: 



2ttM ; 

H 
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with the definitions: 



sin(0-0) = f(9) 



cos{8-4>) - g(0) (7) 
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Fig. 2 - Parametric curve for an infinite cylinder with uniaxial anisotropy. We choose the normalised 
anisotropy constant 2 J^ 2 = 0.5 and several reduced radii S =1,2 and 100. The perfect Stoner- 
Wohlfarth astroid occurs when S = 0. 



g(6) = 2{N a sin 2 9 + N c cos 2 9) - a ^-^(3 cos 2 6-1)- 2 ^ 2 (5 cos 2 6-1) sin 2 6 

2irM s 2ttM s 

(8) 

, 2 

The parameter a = with k = ^ according to Aharoni's [9] work, is a monotonically 
decreasing function of the aspect ratio of the ellipsoid like the parameter q (see Introduction). 
As an example, when S = 2.06 the Nickel nanowire radius is about 41.2 nm (see Figs. 3 and 
4). 

Aharoni's approach is based on an isotropic approximation of the inhomogeneous compo- 
nents of the curling magnetization. Following Muller and Goldstein [16], we have developed 
a different set of curling equations in presence of fourth-order uniaxial that is free from the 
isotropic approximation as well as in presence of cubic anisotropy. In each case we get an 
appropriate g{9) function as given below. 

In our fourth-order uniaxial case [15], g{6) is given by: 

g(9) = 2(N a sin 2 6 + N c cos 2 6) -a - — + 2 ^ cos 29 - 2 ^ 2 2 cos 2 9(3 sin 2 9 - cos 2 9) (9) 

2nM s 2itM s 

whereas in the sixth-order cubic anisotropy [15] case, it is given by: 

g(6) = 2(N a sin 2 6 + N c cos 2 9) - a ^-^(1 - 5 cos 2 9 sin 2 9) sin 2 6 cos 2 6 (10) 

W V ; 2ttM 2V ; 4ttM 2 V ' 
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Equations 8 (along with the different definitions of the g{6) function) can be used to extract 
the components of the switching field h x ,h z as before: 



= 6 — tan 



/(*) 



.9(0). 

h x = sin(%(0) - cos(0)/(0) 
h z = cos(%(0) + s in(0)/(0) 



(11) 



We use the same parameters as Wegrowe et al. [17], K\ = 2. Ox 10 5 erg/cm 3 , K2=0 erg/cm 3 
and compare them to the Nickel bulk values at room temperature. 
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Fig. 3 - Geometric signature in the Wegrowe et alcase [17] (Ki = +2.0 x 10 5 erg/cm 3 whereas 
K2 — 0. erg/cm 3 ) for three aspect ratios m = c/a = 3, 6 and 100. The reduced radius in all cases is 
S = 2.06. 




Fig. 4 - Geometric signature for bulk Ni at room temperature K\ — — 4.5x 10 4 erg/cm 3 , K2 = 2.3 x 10 4 
erg/cm 3 ) for three aspect ratios m = c/a = 3, 6 and 100. The reduced radius is S — 2.06. 
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Fig. 5 - Comparison of least-squares fitting and measurements (+) on all Nickel nanowires in the 
uniaxial case. Geometric signature for the Ni nanowire experimental results with four sets of diam- 
eters: 15, 50, 80 and 100 nm. The fitting values are respectively: M a = 195 emu/cm 3 , ifi=-223420 
erg/cm 3 and K^= 73198 erg/cm 3 for 15 and 50 nm and M a = 56 emu/cm 3 , 7^1=6035 erg/cm 3 and 
if 2 =-7255 erg/cm 3 for 80 and 100 nm diameters. The reduced radius is 5=0.728, 2.427, 3.883 and 
4.854 respectively. The framing curves are obtained from the fitted Ki, if 2 and the bulk value of M s . 



Our experimental results for the four diameter nanowires are fitted with a least-squares 
method to the fourth-order uniaxial case (in the Aharoni and our cases) as well as to the sixth 
order cubic anisotropy case and are displayed in Figs. 5 and 6 for Nickel and the same is done 
for Cobalt in the fourth-order uniaxial case only (see Fig. 7). 

Recalling that for Nickel, the room temperature bulk values are: M 8 — 485 emu/cm 3 , 
cubic [15] anisotropy constants: K\ = —4.5 x 10 4 erg/cm 3 , K2 — 2.3 x 10 4 erg/cm 3 and for 
Cobalt, the room temperature bulk values are: M s — 1400 emu/cm 3 , uniaxial [15] anisotropy 
constants Ki=-9 xlO 5 erg/cm 3 and K 2 ~ -2 xlO 5 erg/cm 3 , the least-squares fitting results 
show that the geometric approach to reversal modes (SW and curling) allows to have a good 
insight regarding the nature of the reversal modes despite the fact some of the fitting values 
depart substantially from the bulk values. 
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Fig. 6 - Comparison of least-squares fitting and measurements (+) on all Nickel nanowires in the 
cubic anisotropy case [15]. Geometric signature for the Ni nanowire experimental results with four 
sets of diameters: 15, 50, 80 and 100 nm. The fitting values are respectively: M s = 195 emu/cm 3 , 
_ftTi=-223420 erg/cm 3 and K 2 = 73198 erg/cm 3 for 15 and 50 nm and M s = 56 emu/cm 3 , ifi=6035 
erg/cm 3 and i\ 2 =-7255 erg/cm 3 for 80 and 100 nm diameters. The reduced radius is 5=0.728, 2.427, 
3.883 and 4.854 respectively. The framing curves are obtained from the fitted K\, K2 and the bulk 
value of M s . 



Discussion and Conclusion. - In conclusion, by means of theoretical studies and experi- 
mental measurements, we have investigated the reversal processes in ferromagnetic nanowires. 
Our systematic studies of the effect of the nanowire radius show that the magnetization re- 
versal mechanism is strongly influenced by its value. Two rotation modes are considered as 
the most important: coherent rotation and curling. Good agreement between the measured 
magnetic properties of Ni nanowires and the theoretical calculations is obtained. However, 
further experimental work remains to be done in order to observe this transition. 

APPENDIX. - The establishment of the geometric representation is based on the fol- 
lowing steps: 
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Fig. 7 - Comparison of least-squares fitting and measurements (+) on all Cobalt nanowires. Geometric 
signature for the Co nanowire experimental results with four sets of diameters: 15, 50, 80 and 100 
nm. The fitting values are: M a = 87 emu/cm 3 , Ki=-11333 erg/cm 3 and K 2 = 11318 erg/cm 3 for 15 
nm, M s = 105 emu/cm 3 , /fi=-11691 erg/cm 3 and K 2 = 7268 erg/cm 3 for 50 nm, M a = 44 emu/cm 3 , 
ill =-9912 erg/cm 3 and K 2 = 22826 erg/cm 3 for 80 nm, and M 3 = 79 emu/cm 3 , K x =-90674 erg/cm 3 
and if 2 =28643 erg/cm 3 for 100 nm diameter. The reduced radius is 5=0.065, 0.262, 0.175 and 0.396 
respectively. The inset in the 100 nm case shows the detailed quality of the fit. The framing curves 
are obtained from the fitted K\ , K 2 and the bulk value of M s . 



1. The equilibrium magnetization in the uniform case is first established. This entails 
writing the energy as E = Ek + E m + Eh where Ek is the anisotropy energy, E m the 
magnetostatic (demagnetization) energy and Eh the Zccman energy. The derivative 
75§ = gives the equilibrium uniform magnetization Mq from the energy minimization. 

2. The perturbation of the equlibrium (uniform) magnetization Mo by a non- uniform vari- 
ation em(r) leads to the space dependent total magnetization M(r) = M + em(r) 
(see fig. 1) with the condition of local orthogonality em(r) • M = 0. From the non- 
uniform magnetization, M(r) one calculates the average value of the energy E = 
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Ek + E m + Eh + Ex over the nanowire where one explicitly introduces the exchange 
energy Ex stemming from non-uniformity. The total energy is now an even function of 
e from symmetry consideration and is written as: E(9, e 2 ). The minimization of E(9, e 2 ) 
with respect to the perturbing amplitude e 2 i.e. |^ = gives the second set of equations 
that leads to the sought parametric representation. 

For an Infinite cylinder, Shtrikman and Treves [18] were the first to derive the solution of 
the non-uniform magnetization in the (u, w, y) system (see fig. 1) as given by: 
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with A 

first positive maximum of J\{x) is Xi 
From the series representation of the Bessel function Ji(Xx): 



the positive roots of the equation dJ ^^ = 0. The first solution (n = 1), that is the 

1.841 (see for instance Abramowitz and Stegun [10]). 



we get: 
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hence the Ishii [13] approximation for the Ji(Aix) function: 
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that he used to evaluate the average energy over the cylinder. 

The magnetization components to be considered are Mi which are in the rotated frame 
(u,w,y) (see fig. 1) with u along the equlibrium magnetization Mq. The magnetostatic 
energy 2wNijMiMj ,i, j = 1,2,3 averaged over the volume of an infinite cylinder of radius 
a is given approximately by the average of 2n NijMi(r)Mj(r) [19]. This means taking the 
demagnetization coefficients of the full cylinder and letting the magnetization M(r) carry 
the spatial dependence. The averaging in the (u,w,y) system yields: 
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The average exchange energy is given by: 
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The average anisotropy and Zeeman energies are evaluated likewise such that: 
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